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Abstract 

Using the recently derived representation for the polarization tensor in (2+l)-dimensional space- 
time allowing an analytic continuation to the entire plane of complex frequencies, we obtain simple 
analytic expressions for the reflection coefficients of graphene at nonzero Matsubara frequencies. In 
the framework of the Lifshitz theory, these coefficients are shown to lead to nearly exact results for 
the Casimir free energy and pressure between two graphene sheets. The constituent parts of large 
thermal effect, arising in the Casimir interaction at short separations due to an explicit parametric 
dependence of the polarization tensor on the temperature and an implicit dependence through a 
summation over the Matsubara frequencies, are calculated. It is demonstrated that an explicit 
thermal effect exceeds an implicit one at shorter separations, both effects are similar in magnitudes 
at moderate separations, and that an implicit effect becomes a larger one with further increase of 
separation. Possible applications of the developed formalism, other than the Casimir effect, are 
discussed. 

PACS numbers: 78.67.Wj, 42.50.Lc, 65.80.Ck, 12.20.-m 
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I. INTRODUCTION 


Graphene has already assumed great importance as a material possessing unusual me 


Q- 


chanical, electrical and optical properties of high promise for various applications HI. These 
properties are due to 2D character of graphene. In contrast to standard 3D materials, quasi¬ 
particles in pristine (undoped) graphene are massless Dirac fermions which possess the linear 
dispersion relation but move with the Fermi velocity rather than with the velocity of light 
lj, 2], Taking into account that graphene is considered as a potential element of rnicro- 
and nanoelectromechanical devices, much attention is attracted to forces acting between 
two graphene sheets and between a graphene sheet and some element made of an ordinary 
material (metallic, dielectric or semiconductor). These forces are caused by the zero-point 


and thermal fluctuations of the electromagnetic field and are callec 


the van der Waals or 


.a. 


The fluctuation 


Casimir forces depending on the range of separations considered 
induced forces are known also under the generic name of dispersion forces. 

The van der Waals and Casimir forces acting between two graphene sheets and between 


a graphene sheet and a 3D material were studied theoretical 


ly by many authors in the 


15]). The most important 


framework of various approaches (see, for instance, Refs, 
breakthrough was made in Ref. [8], where it was shown that for graphene the thermal 
effect in the Casimir force becomes crucial at by an order of magnitude shorter separations 
than for ordinary materials. This result was obtained using the longitudinal density-density 
correlation function in the random-phase approximation, and the relativistic effects were 
shown to be not essential. 

The most straightforward description of the van der Waals and Casimir forces in layered 
systems including graphene is given by the Lifshitz theory j^, 4, 16] with reflection coefficients 
expressed via the polarization tensor of graphene in (2+l)-dimensional space-time. The 
explicit analytic expressions for the components of this tensor were obtained in Refs. HQ 


at zero and nonzero temperature, respectively. Using the values of the polarization tensor 
of graphene at the imaginary Matsubara frequencies, a lot of numerical computations of 
the Casimir and Casimir-Polder forces between graphene and different materials and atoms 
has been performed 


17 


-271 ]. The formalisms exploiting the polarization tensor and the 


density-density correlation function have been compared 


28). It was shown 28] that these 


formalisms are equivalent if one uses the exact temperature-dependent longitudinal and 
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transverse density-density correlation functions. Computations using the polarization tensor 
confirmed the existence of large thermal effect in the Casimir interaction of graphene sheets 

8]. Furthermore, the results computed using the 


at short separations, as predicted in Ref. 
polarization tensor of graphene at nonzero temperature 


M 


30| were found to be in very good 


agreement with first measurement of the gradient of the Casimir force between a Au-coated 


sphere and a graphene-coated substrate by means of a dynamic atomic force microscope 




The computational results obtained using the hydrodynamic model of graphene la la 


11 


were shown to be excluded by the measurement data over a wide range of separations 


31J. 


26] 


32]. 


In this paper, we investigate an origin of the large thermal effect in the interaction of two 
graphene sheets at short separations from a few nanometers to a few hundred nanometers. 
Our aim is to determine the relative contributions to the Casimir free energy and pressure 
due to an explicit dependence of the polarization tensor on the temperature as a parameter 
and an implicit dependence on the temperature through a summation over the Matsubara 
frequencies. For this purpose, we circumvent cumbersome computations using the exact 
expressions for the polarization tensor and find a transparent way of performing simple 
analytic calculations in powers of a natural small parameter as far as possible. It is significant 
that the polarization tensor of graphene at nonzero temperature obtained in Ref. 18] is not 
well adapted for analytic calculations. The point is that it is valid only at the discrete 
imaginary Matsubara frequencies and its immediate continuation to the entire imaginary 
frequency axis does not satisfy necessary physical conditions 33]. Another form of the 
polarization tensor of graphene, valid over the whole plane of complex frequency, including 
the imaginary and real frequency axes, was found in Ref. ( 33 1. At the imaginary Matsubara 
frequencies the polarization tensor of Ref. 33] takes the same values as the one of Ref. 18], 
but presents better opportunities for analytic calculations due to much simpler mathematical 
structure. 


Here, we use the polarization tensor in the representation of Ref. 33] and obtain simple 


asymptotic expressions for both the thermal correction to this tensor and for the reflec¬ 
tion coefficients of the electromagnetic oscillations from a graphene sheet at the Matsubara 
frequencies. The developed formalism is applied to calculate the Casimir free energy and 
pressure between two graphene sheets at room temperature. We compare our results with 


previously made 


22] numerical computations using exact expressions for the polarization 


tensor in the representation of Ref. 18]. We demonstrate that the error of our asymptotic 
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results is only a small fraction of a percent over the whole region of physically meaningful 
values of parameters. Then we investigate the origin of large thermal effect at short sep¬ 
arations. It is shown that the explicit and implicit contributions to the dependencies of 
the Casimir free energy and pressure of graphene on the temperature are equally important 
parts of the large thermal effect. In so doing, the major contribution to the thermal effect 
due to an explicit temperature dependence of the polarization tensor originates from the 
zero-frequency Matsubara term. 

The structure of the paper is the following. In Sec. II we introduce the main notations, 
formulate our perturbation approach and obtain simple asymptotic expressions for the po¬ 
larization tensor of graphene and reflection coefficients. Section III is devoted to the inves¬ 
tigation of large thermal effect in the Casimir free energy of two parallel graphene sheets. 
In Sec. IV the same is done for the Casimir pressure. Section V contains our conclusions 
and discussion. 


II. SIMPLE ASYMPTOTIC EXPRESSIONS FOR THE POLARIZATION TEN¬ 
SOR OF GRAPHENE AND REFLECTION COEFFICIENTS 


In the framework of the Dirac model, the reflection coefficients of the electromagnetic os¬ 
cillations from graphene can be expressed via the polarization tensor of graphene in (2+1)- 
dimensional space-time. To calculate the Casimir free energy and pressure between two 
graphene sheets one needs the reflection coefficients calculated at the pure imaginary Mat¬ 


subara frequencies 


3 ,y, 


22 ] 


rTM k_ l) = 
tte^i , k_ l) = 


qiU 00 (i^i,k A 


qiU 00 (i^i, k ± ) + 2 hkY 

n(?6 ,k±) 


(i) 


n(i&, k_ l) + 2 hk\qi 

for two independent polarizations of the electromagnetic held, transverse magnetic (TM) 
and transverse electric (TE). Here, = 27 t ksTl/h, where ks is the Boltzmann constant, T 
is the temperature and / = 0, 1,2, ..., are the Matsubara frequencies, k± is the magnitude 
of the projection of the photon wave vector on the plane of graphene, qi = (k\ + ^f/c 2 ) 1 / 2 , 
n mn with m, n = 0, 1, 2 are the components of the polarization tensor of graphene in 
(2+l)-dimensional space-time, and the quantity n is defined as 


n(?6, k_ l) = k\ n tr (i&, k_ l) - q?Tloo{i£,i, k A 


( 2 ) 
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where n tr = II//. Note that the reflection coefficients (JT]) can be equivalently rewritten via 
the longitudinal and transverse density-density correlation functions, or conductivities, or 
polarizabilities or nonlocal dielectric permittivities of graphene 28]. 

The components of the polarization tensor of graphene at temperature T can be presented 
in the form 

n 0 o(h£z, k_ l) = II w(i£i, k ±) + A T noo(i&, k _ L ), 

n(?6, k ± ) = n {0) (i6, k ± ) + A r n(^, k ± ), (3) 

where n^, iR 0 ) are the respective quantities calculated at zero temperature, but with 
continuous frequencies £ replaced by the Matsubara frequencies £;, and A r IIoo, A^II are the 
thermal corrections to them. 

For a pristine (gapless) graphene under consideration in this paper, the components of 
the polarization tensor at T = 0 K take especially simple form 17] 


nS ) « I ,U) = 


irahk 2 

a 


II (0) (i£,fci) = 7 rahk 2 ± q h 


where a = e 2 /(hc) ~ 1/137 in the fine-structure constant, 


(4) 


(5) 


and vp ~ c/300 is the Fermi velocity. The respective reflection coefficients (JT]) are given by 




iraqi 


4e(^ 6, fc±) = 


iraqi 


( 6 ) 


7rag; + 2g, Iraqi + 2q t 

For the thermal corrections to the polarization tensor we use the expressions obtained in 
Ref. 33] 

3 dq ± 


1 Qafic 2 

A r n 00 &t,k ± ) = -^F_ 


e hcqj_/(k B T) _|_ -[_ 


X 


1 - —ReW^i,k ± ,q ± ) 
Qi 


(7) 


, , 16 ahc 2 

ArII(«£i, k±) = -2— 


x i + ^ 


dg± 

e hcq ± /(k B T) _|_ l 


^-Re 


W^,k ± ,q ± ) 


Re VF(6, k±, q±) — 

Here, the integration is performed over the magnitude of a two-dimensional wave vector of 
the electronic excitation and the quantity W is defined as 


W (£z, k±, q±) = qf-4qi + 




( 8 ) 
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The expressions in Eq. (ED have the same values as the respective more complicated expres¬ 


sions in Ref. 


18], but, in contrast to the latter, can be immediately analytically continued 


to the whole plane of complex frequency. 

It is possible to simplify Eq. (O taking into account that 


Re IE (6, k±, q±) = 


Re 


1 


v/2 


(qf - 4q\) 2 + 16 ^ 2 g± + qf - 4 q]_ 


1/2 


ReIE(6, k±,q±) 


W(^k ± ,q ± ) 


(9) 


Using Eqs. (J5]) and ([9]) and introducing the new integration variable u = 2 q±/qi, we rewrite 
Eq. © in the form 

8 ahc 2 qi f°° du 


A T n 00 (i6, k ± ) — 


e B ‘ u + 1 


x \ 1 — 


72 


(1 + MT _ 4 !M/i + 1 _„ 2 


c 2 qf 


1 / 2 ' 


( 10 ) 


A TT ft 1 \ 8ahc2 U 

A T n(?6, k±) = -2— 


du 


e B ‘ u + 1 


x 


x 


c 2 + V 2 


(1 + u 2 ) 2 — 4 


k 2 , u 2 


J F n '± 


c 2 qf 


+ 1 -u 2 


1/2 


1 - 


v\k\ 


c 2 qf {1 + u 2 ) 2 — 4 


Vp “ 2 

c2 4f 


where Bi = hcqi/(2k B T). 

Now we consider approximate expressions for the quantities (flOD applicable at any l > 1 
for calculations of the Casimir force using the Lifshitz formula. For this purpose we note that 
at room temperature T = 300 K the first Matsubara frequency £i ~ 2.4 x 10 14 rad/s. Taking 
into account that the characteristic photon wave vector giving the major contribution into 
the Casimir free energy and pressure is equal to k± = 1/(2a), where a is the separation dis¬ 
tance between the plates [l, 34], one concludes that the natural parameter entering Eq. (fTOD 
can be considered as a small one at all separations a > 10 nm, 


4 v 2 F k\ 
c 2 qf 


< 


4 v 2 k 2 1 

c 2 q( 


( 11 ) 


Thus, at a = 10 nm this parameter is less than 0.17, at a = 20 nm less than 0.04 and further 
decreases with increase of separation. 
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Now we expand the square roots in Eq. (TTOT) in powers of the small parameter (HIT) and 
preserve only the zero and first order contributions. Furthermore, from Eq. m we have 
Vpk“]_ -C £1 and, thus, Bi ~ h^i/{2k B T) = 7 tI. Then from Eq. (HH one obtains 


ryfik 2 

A T n 00 (i£ h k ± ) = 


A r II (i£i,k±) = ahkpqiYi, 


( 12 ) 


where we have introduced the notation 


Y, = 4 


du u~ 


e nlu _|_ l 1 -|- 


(13) 


By combining Eqs. (j4j) and (TT2ll in accordance with Eq. (J3]) , we arrive to the following 
simple asymptotic expressions for the polarization tensor of graphene 


„ . „ , . ahk 2 , . 

n 0 o«/, k±) = —+ Yi), 
Qi 

n(i&, k ± ) = ahk]_q l ('K + Y t ). 


(14) 


Substituting Eq. (T14T) in Eq. (JT]), one obtains simple equations for the reflection coefficients 

aq^n + Yi) 


^tm(*6i k±) = 
rTs{iCi, k ±) = 


aqi(n + Y t ) + 2 qi 

aqi(n + Yi ) 


(15) 


aqi(n + Yi) + 2 qi 

As can be seen in Eqs. (fT4]l and f[T5]h in both the polarization tensor and reflection coefficients 
the relative contribution of each term with l > 1 is characterized by the ratio Yj/ 7 r. In this 
connection it is instructive that for l — 1, 2, 3, 4, and 5 the quantity Yi/tc is equal to 0.041, 
0.0074, 0.0024, 0.0011, and 5.5 x 10” 4 , respectively, and further decreases with increasing l. 

Recall that Eqs. o and 03 are valid under a condition l > 1. For l = 0 the thermal 
correction to IIoo can be obtained from the first formula in Eq. (TTOl) 


A r n 00 (0, k ± ) = ^(k B T In 2 hVpk ' 


du 


e B ° u + 1 


vT^ 


U‘ 


(16) 


where Bq = hvpk±/(2k B T). In the derivation of Eq. (fT6l) we have taken into account that 


7 T 


a 2 ) 2 = 


u~ 


u < 1, 


u 


1, u > 1. 


(17) 
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(18) 


Using the first formula in Eq. (H|) . from Eqs. ©, © and (TT6|) we obtain 

ac {WksT In 2 + hk±v f[k — 8A(fcj)]} 


r tm ( 0 , k±) — 


2 ’ 


ac{16fc j BTln2 + hk±VF[n — 8X(Vl)]} + 2 hk±Vp 


where 


X(k ± ) = 


du 


-VI 


u 


J 0 e B ° u + 1 

Equation (jT8|) is an exact one. It can be rewritten in an equivalent form 

2hk±Vp 


^tm(0, k±) = 1 - 


(19) 


( 20 ) 


ac {16VbT1ii 2 + hk±vp[TT — 8A(Aq_)]} + 2 hk±Vp 
From Eq. ()20)) it is seen that at zero Matsubara frequency the values of the TM reflection 
coefficient are rather close to unity and 


lim r TM (0, kj_) = 1. 

/ c _|_ —^0 


( 21 ) 


In a similar way, for l = 0 the thermal correction to II can be found from the second 
formula in Eq. (fTOl) 

ATmk±) = _^^i Z(k± ), (22) 

c 

where 


Z(k ± ) = 


du 


u 


Jo e s «“ + lvT^V' 

Using the second formula in Eq. ((4]), from Eqs. ([!]), (j3j) and (T22|i one obtains 

avp[ 7r — 8Z(k±)\ 


^te(0, k_ l) — 


avp[ 7f — 8Z(k±)\ + 2c 


It is easily seen that 


|?"te(0, k± | < a ^ F & 3.8 x 10 5 


(23) 


(24) 


(25) 


and, thus, the TE mode does not contribute into the classical limit 


23 


24] 


III. INVESTIGATION OF THE THERMAL EFFECT FOR THE CASIMIR FREE 
ENERGY 


The Casimir free energy per unit area of two parallel graphene sheets separated by a 
distance a at temperature T in thermal equilibrium with an environment is given by the 


Lifshitz formula 


BE, 


16] 


h T °° / f°° 

T(a,T) = ^~Yl / k ± dk ± {\n[l-r 2 TM (i^,k ± )e- 2 ^] 


1=0 


ln[l-r£ E (i&,k±)e 2aqi }} , 


(26) 
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where the prime on the summation sign indicates that the term with l — 0 is divided by two 
and the reflection coefficients are given in Eq. ©■ 

For convenience in numerical computations we rewrite Eq. (127)1) in terms of dimensionless 
variables 

V = 2q t a, Ci = —. (27) 

c 

Then the Casimir free energy takes the form 

T ) = / ydy (M 1 " r ™«i» y)e~ y ] 

i=o 

+ ln[l — r^ E (iCi, y)e~ y ]} , (28) 


where the reflection coefficients are 


rnssiKhV) 


^te (i£i,y) 


ytl m (iC h y) 

y^oo{Khy) + ~ Ci)' 

_ n (iCi,y) 

fl(iCi,y) + 2y(y*-C?) 


and the dimensionless polarization tensor is given by 


TT 2 a 

‘' rv n t 11 mn ? 

n 

o„3 

n = —n = (y 2 - c, 2 )n„ - /n„„. 


(29) 


(30) 


In terms of the variables (1271) . the dimensionless polarization tensor at zero temperature 
OH) is equal to 


n 0 o(iCi,3/) 
n (iCi,y) = 


na(y 2 - Cf) 

fi 

na{y 2 - Ci)fh 


(31) 


where 

fi = fi{y) = [vW + (1 - v 2 f)Ci] 1/2 (32) 

and Vp = Vf/c is the dimensionless Fermi velocity. The respective reflection coefficients at 
T = 0 K (where continuous frequency is replaced with discrete Matsubara frequencies) are 
obtained from Eq. (1231) 


4m (Ki,y) 
4e (iCi,y) 


iray 

nay +2 ft ’ 
nafi + 2 y' 


( 33 ) 
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The dimensionless thermal correction to the polarization tensor (13111 is easily obtained 
from Eq. (fTOjl 

8 aft r 


Arn 0 o(*0; y ) — 


u 


■ F JO 


e B i u + i 


x < 1 - 


A T n (iCi,y) = 


V2 

8 aft 


' (1+U 2 )2 _ 4 flVziLhi + w 

30 du 


1 1/2 ' 


v“p Jo e B ‘ u + 1 


-C 2 + 


ft 

V2 


x 


(1 + u 2 ) 2 — 4 


W-tf) 


W - O 2 


f 


+ 1 -u 2 


1/2 


vl{y 2 -(?)u 2 


(34) 


where Bi = rcfi/r and r = 47 t ak F T j{hc) is dimensionless temperature. 

Using Eqs. urn and (13ip . we find approximate expressions for the polarization tensor in 
terms of dimensionless variables which are valid at all nonzero Matsubara frequencies 


fiooKi.a) = ° :(i<2 , C,2) (T + y,). 
n 

fl(i(i,y) = a(y 2 -(?)f l ('K + Y l ), 


(35) 


where Y] is given by Eq. (fT3]l . The respective approximate reflection coefficients (fT5ll valid 
at l > 1 take the form 


rTM(i(i,y) 


rrE(i(i,y ) 


ay(7T + Yi) 
ay(n + Y t ) + 2 fj 

_ afi(n + Y t ) 

afi(n + Yi) + 2y' 


(36) 


Finally, the exact reflection coefficients (120P and (1241) valid at zero Matsubara frequency 
in terms of dimensionless variables are given by 


ftm(0 ,y) 

t-te(0,2/) 


2irv 2 F y 

a {8rln2 + nvFy[^ — 8X(?/)]} + 2nv F y' 
av F [n - 8Z(y)\ 
av F [ 7r — 8 Z(y)\ + 2’ 


(37) 


where the quantities X(y) and Z{y) are defined in Eqs. (fT9l) and (j23j) taking into account 
that Bo = nv F y/ t. 
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Now we are in a position to determine the accuracy of our approximate expressions and 
to investigate the origin of large thermal effect in the Casimir free energy between two 
graphene sheets at short separations. For this purpose we calculate the Casimir free energy 
(j28|) per unit area of two graphene sheets in three different ways. First, we compute the 
exact free energy T ) using Eqs. (125|) and Q29J) , where the polarization tensor is given by 
the sum of respective expressions in Eqs. (l3Tj) and (1341) . Then, we perform calculations of the 
approximate free energy ^^(a^T) by using the exact reflection coefficients (1371) at Co = 0 
and our asymptotic reflection coefficients (1361) at all Q with l > 1. Finally, we calculate 
the free energy J r<2 ^ ) (a, T) using the approach of Ref. 18j, i.e., by using the exact reflection 
coefficients (1HTD at Co = 0 and coefficients (fool) at all l > 1. Note that the approximate 
free energy T (l \a, T ) takes an exact account of the explicit temperature dependence of the 
polarization tensor in the zero-frequency term of the Lifshitz formula. In all terms with 
l > 1 the explicit temperature dependence is taken into account approximately using our 
asymptotic approach. The approximate free energy J rl ^ 2 \a,T) takes into account an explicit 
temperature dependence only in the zero-frequency term and discards it in all terms with 
l > 1 . 


To compare the accuracies of both approximate calculation approaches, in Fig. QJa) we 
plot the quantity 

S ^ T) = T{a ’ T) -^ (a ’ T) (») 
7 ( 0 , 7 ) 

at T — 300 K as a function of separation by the lines 1 and 2 in the region from a = 5 to 
250 nm (k = 1, 2) for the first and second approaches, respectively. As is seen in Fig. [D(a) 
(line 1), the error of our asymptotic approach is negative and the maximum of its magnitude 
equal to 0.17% is achieved at the shortest separation a = 5 nm. With increase of separation to 
10 and 30 nm the error magnitude decreases to 0.12% and 0.02%, respectively. At separations 
a > 35 nm our calculation approach becomes practically exact. 

Note that the characteristic value of y giving the major contribution to the Casimir 
free energy (l28j) is y — 1. Taking this into account, the inequality (ITT]) reduces to r > vf- 
This inequality determines the region of separations and temperatures where our asymptotic 
approach is applicable for calculations of the Casimir free energy and pressure (for the latter 
see Sec. IV) in the framework of the Lifshitz theory. According to the results of Sec. II, at 
T = 300 K it should work well at a > 10 nm, where the parameter (TTIT) is less than 0.17. 
From Fig. [D(a) (line 1) it is seen that the proposed asymptotic approach leads to rather 


11 




accurate result even at a = 5nm where the parameter (ITT]) is equal to 0.6. One should take 
into account, however, that the Dirac model of graphene is only applicable at energies below 
a few eV. This questions the possibility of using it for theoretical description of the Casimir 


force at separations below approximately 20 or even 30 nm 


32]. 


The second approximate way of calculations turned out to be less accurate than the 
first one. As is seen from line 2 in Fig. QJa), in this case the calculation error is positive. 
It is equal to 0.36% at a = 5 nm, achieves the largest value of 0.78% at a = 20 nm and 
then decreases to 0.49%, 0.21%, and 0.05% at separations 50, 100, and 250 nm, respectively. 
Taking into account, however, that in the second way of calculations the error remains to be 
below 1%, we conclude that the major contribution of an explicit temperature dependence 
of the polarization tensor to the Casimir effect originates from the zero-frequency Matsubara 
term. 


Now we discuss an origin and sources of 
arising for graphene at short separations 


die large 


18 


20 


hermal effect in the Casimir free energy 
E|. For this purpose in Fig. [2]a) we 
plot the magnitudes of the Casimir free energy normalized to the quantity A = ksT /(87ra 2 ) 
at T = 300 K [see Eq. (128|) ] as functions of separation in the region from 5 to 250 nm. The 
upper line (J 7 ) was computed at T = 300 K using the Lifshitz formula fl28jl and the reflection 
coefficients (l29jl with the exact polarization tensor defined as a sum of respective expressions 
in Eqs. (13T1) and (134|) . The nearly coinciding line [T = J-A)] j s obtained also by using our 
asymptotic reflection coefficients (j36l) at l > 1 and (l37l) at l = 0. The intermediate line [J 7 ^] 
was computed at T = 300 K using the Lifshitz formula (1281) and the reflection coefficients 
(l33lh i.e., using the polarization tensor at T = 0, where continuous frequency was replaced by 
the discrete Matsubara frequencies. In this case an explicit dependence of the polarization 
tensor on T was discarded. Finally, the lowest line (T'o) presents the Casimir free energy 
at T = 0 (i.e., the Casimir energy E(a) per unit area of two graphene sheets) which, by 
definition, does not contain any thermal effect, either explicit or implicit. 

As can be seen in Fig. [2](a), there is large implicit thermal effect for two graphene sheets 
at short separations due to a summation over the discrete frequencies. The magnitude of this 
effect is characterized by the difference between the intermediate and lowest lines [for better 
visualization in Fig. [2](b) the region from 5 to 30 nm is shown separately on an enlarged 
scale]. The difference between the upper and intermediate lines characterizes the magnitude 
of the explicit thermal effect arising from the dependence of the polarization tensor on the 
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temperature as a parameter [see also Fig. [2^b) at shorter separations]. As to the difference 
between the upper and lowest lines, it represents the total thermal correction to the Casimir 
energy. As is seen in Fig. [21 at moderate separations the explicit and implicit contributions 
to the thermal effect are nearly equal. 

In Table I we illustrate the quantitative relationships between different constituent parts 
of the thermal correction and the role of the total thermal correction in the Casimir free 
energy at T — 300K. In columns 2, 3, and 4 the fractions of explicit, ([T 7 ! — |J r(0 ^|)/|J r |, 
implicit, — |J 7 0 |)/| J r|, and total ([J 7 ! — |^ r o|)/|-? 7 |, thermal corrections are presented 

at different separations indicated in column 1. The quantities \T\, IJ 7 ^ 0 ^ and |J 7 0 | = \E\ 
are those plotted by the upper, intermediate and lowest lines in Fig. [21 respectively. As is 
seen in Table I (column 4), the thermal effect in graphene quickly increases with increase of 
separation and provides more than 80% contribution to the magnitude of the Casimir free 
energy at separations above 100 nm. In so doing at separations below 200 nm the TM mode 
contributes more than 99% of the free energy and more than 99.9% at a > 200 nm. An 
explicit thermal effect (column 2) contributes more than one half of the thermal correction 
at separations below 155 nm and at larger separations becomes less than an implicit thermal 
effect. However, our calculations show, that the magnitudes of both effects are similar and, 
thus, both of them must be taken into account in computations. We note also the role of 
the Matsubara term with 1 = 0. At separation distances of 30, 100, 200, and 1000 nm it 
contributes 69.3%, 94.4%, 98.4%, and 99.94% of the Casimir free energy, respectively. 

As was mentioned in Sec. I, the thermal effect in the Casimir force for graphene becomes 
dominant at much shorter separations than for plates made of metallic or dielectric materials. 
For example, for ideal metal plates the relative thermal correction to the Casimir free energy 
at T = 300 K is equal to ([J 7 ! — |.7 7 o|)/|.7 7 | = 3 x 10 -5 and 2.6 x 10~ 2 at separation distances 
a = 100 and 1000 nm, respectively. This should be compared with Table I, where the 
respective relative thermal effects for graphene are indicated as 0.7922 and 0.9788, i.e., by 
the factors of 2.6 x 10 5 and 37.6 larger. 

The physical reason for such a big difference between graphene and other materials is 
the following. For ordinary metals and dielectrics the thermal effect becomes dominant 
when the contribution of all Matsubara terms with l > 1 becomes exponentially small 
due to the exponential factors in Eq. (1261) .At T = 300 K this happens at separations 
above 5 or 6 /jrn. In doing so, the remaining zero-frequency Matsubara term alone describes 
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the classical thermal Casimir effect. As to graphene, at short separations above 100 nm, 
where the exponential factors in Eq. (126]) are not small yet, the reflection coefficients in 
all Matsubara terms with l > 1 become very small making dominant the classical thermal 
contribution of the zero-frequency term. 


IV. INVESTIGATION OF THE THERMAL EFFECT FOR THE CASIMIR PRES¬ 
SURE 


The Lifshitz formula for the Casimir pressure between two parallel graphene sheets is 
given by 0,0, Q 


P(a ,T) = --^V / qik ± dk ± {[r T l l (i£i,k ± )e 2aqi - 1] 1 

vr ' ,/n 


1=0 


+ [ r TEKi.U)e 2 “® - 1] ’}, 


(39) 


where the exact reflection coefficients are written in Eq. ([!]) with the polarization tensor 
presented in Eqs. ([3]), (J4J) and dZ|). 

In terms of the dimensionless variables (1271) convenient in numerical computations 
Eq. (139|) takes the form 


kuT 


p(a,T) = --5-J / y 2 dy { I'tmUO, y)e y - 1] 1 


8rf — •'£, 

+ [ril«i,y)e v - l] -1 }. 


(40) 


where the exact reflection coefficients can be found in Eq. (1291) and the dimensionless polar¬ 
ization tensor is presented in Eqs. (13TT) and (1341) . The approximate reflection coefficients at all 
nonzero Matsubara frequencies (Z > 1) using our asymptotic approach are given in Eq. (136|) . 
At zero Matsubara frequency the exact reflection coefficients are written in Eq. (l37l) . 

First, we determine the accuracy of our approximate method in application to calculations 
of the Casimir pressure. To do this, we calculate the Casimir pressure (TIU1) in three different 
ways already discussed in Sec. III. As a result, we obtain the exact Casimir pressure P(a, T), 
the approximate Casimir pressure P <Kl \a } T) computed using our asy mptotic approach and 
the approximate pressure P ly2 \a 1 T ) using the approach of Ref. 18] discarding an explicit 
dependence of the polarization tensor on T in Matsubara terms with l > 1. The accuracies 
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of both approximate methods are characterized by the quantities 


AP (fc) (a,T) 


P(a,T) - P( k \a,T) 
P(a,T) 


( 41 ) 


In Fig. [T](b) we plot the quantities APd) an d §p( 2 ) at T = 300 K as a function of separation 
by the lines 1 and 2, respectively. As is seen in Fig. [T|(b) (line 1), the error of our asymptotic 
approach in calculations of the Casimir pressure is negative and the maximum of its mag¬ 
nitude is equal to 0.18% at a = 5nm. At separations of 10 and 30 nm the error magnitude 
decreases to 0.16% and 0.01%, respectively, and becomes nearly zero at a > 35 nm. 

For the second approximate way of calculations [the line 2 in Fig. mb)] the error is positive. 
It is equal to 0.22% at a = 5 nm, achieves the largest value of 0.87% at a = 30 nm and then 
decreases to 0.69%, 0.35%, and 0.09% at separations 50, 100, and 250 nm, respectively. Thus, 
our asymptotic expressions for the polarization tensor of graphene and reflection coefficients 
lead to more exact results than the approximate method discarding an explicit temperature 
dependence of the polarization tensor at all l > 1. Nevertheless, the main contribution to 
the explicit thermal effect in the Casimir pressure, similar to the free energy, originates from 
the zero-frequency Matsubara term. 

To investigate an origin and sources of the large thermal effect, in Fig. |3](a) we plot 
the magnitudes of the Casimir pressure normalized to the quantity B = ksT/^Sna 3 ) at 
T = 300 K [see Eq. (128]) ] as functions of separation. The upper, intermediate and lowest 
lines in Fig. [3](a,b) were obtained in the same ways as the respective lines in Fig. [2](a,t)) , i.e., 
the upper line [P = P^ 1 )] was computed at T = 300 K exactly by using Eqs. (l29lh (T3TT) . (l34j) 
and (]40lh the intermediate line [P^] at T — 300K by Eqs. (]33l) and (]40lh and the lowest 
line (Po) is the Casimir pressure computed at T — 0. 

From Fig. |3](a) it is seen that, similar to the free energy, there is large implicit thermal 
effect in the Casimir pressure between two graphene sheets at short separations, character¬ 
ized by the difference between the intermediate and lowest lines. [The separation region 
from 5 to 30 nm is shown on an enlarged scale in Fig. |3](b) .] The explicit thermal effect 
in the Casimir pressure is shown in Fig. |3](a,b) by the difference between the upper and 
intermediate lines. This contribution to the total thermal correction (the latter is shown by 
the difference between the upper and lowest lines) originates from the dependence of the 
polarization tensor on the temperature as a parameter. Again, at moderate separations the 
explicit and implicit thermal effects contribute to the total thermal correction nearly equal. 
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The quantitative results are listed in Table II. The first column of this table contains 
the values of separation and the columns 2, 3, and 4 the values of (|P| — |P (0 )|)/|P|, 
(|p(°)| _ \P 0 \)/\P\, and (|P| — |P 0 |)/|P|, where the magnitudes of the Casimir pressure |P|, 
|Pf°)|, and |P 0 | are the exact ones (the upper line), the ones computed at T — 300K using 
the zero-temperature polarization tensor (the intermediate line), and the ones computed at 
T — 0 (the lowest line), respectively. Similar to the case of the free energy, the thermal effect 
in the Casimir pressure quickly increases with increase of separation (column 4). At sepa¬ 
rations above 150 nm the thermal effect contributes more than 80% of the Casimir pressure 
between two graphene sheets. The contribution of an explicit thermal effect in the Casimir 
pressure (column 2) is larger than that of an implicit at all separations below or equal 
325 nm. At separation distances 30, 100, 200, and 1000 nm the zero-frequency Matsubara 
term contributes 55.7%, 89.9%, 96.9%, and 99.86% of the Casimir pressure, respectively. 
Similar to the case of the Casimir free energy, the thermal effect for the Casimir pressure 
becomes dominant at much shorter separations than for metallic and dielectric materials. 
Thus, for two ideal metal plates at separation distances a = 100 and 1000 nm one has 
(|P| — |P 0 |)/|P| = 1.57 x 1CT' and 1.57 x 1CT 3 , respectively. From Table II we conclude 
that for graphene the relative thermal effect at the same respective separations is larger by 
the factors of 5.6 x 10' and 6.2 x 10 2 . As to the relative role of the TM and TE modes in 
the total Casimir pressure, the TM mode remains dominant and provides more than 99% 
of the pressure magnitude. In precise computations, however, both modes should be taken 
into account. 


V. CONCLUSIONS AND DISCUSSION 


In the foregoing we have obtained simple asymptotic expressions for the polarization 
tensor of graphene and reflection coefficients at nonzero Matsubara frequencies suitable for 
calculation of the Casimir free energy and pressure between two parallel graphene_sheets. 
This has been made possible by applying the new form of the polarization tensor 33] valid 
not only at the imaginary Matsubara frequencies, but over the entire plane of complex 
frequency. In this paper, the analytic asymptotic expressions are derived by means of per¬ 
turbation expressions in powers of a natural parameter which becomes small over the wide 
range of physically interesting separations and temperatures. The respective expressions for 
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the Casimir free energy and pressure between two graphene sheets take into account both 
an explicit dependence of the polarization tensor (or, equivalently, conductivity, or dielectric 
permittivity, or density-density correlation function) on the temperature as a parameter and 
an implicit temperature dependence through a summation over the Matsubara frequencies. 

Using the obtained asymptotic expressions for the reflection coefficients, we have com¬ 
pared the Casimir free energy and pressure calculated in the framework of our formalism 


22] calculated numerically using the polarization tensor of 


with the previously known ones 

Ref. Q. ft was shown that our asymptotic approach leads to very accurate results (the 
maximum errors are equal to 0.17% and 0.18% at the separation of 5nm for the Casimir 
free energy and pressure, respectively). These errors decrease with increase of separation. 
At all separations above 35 nm the asymptotic Casimir free energies and pressures coincide 
with the exact ones. 

We have investigated an origin of large thermal effect arising in the Casimir interaction 
between graphene sheets at short separations [sj]. We have shown that the total thermal 
corrections to the Casimir energy and pressure at zero temperature consist of two parts. 
The explicit part is determined by the parametric dependence of the polarization tensor on 
T, whereas the implicit one arises through a summation over the Matsubara frequencies. 
According to our results, in both the Casimir free energy and pressure, an explicit thermal 
effect exceeds an implicit one at shorter separations. At moderate separations both parts of 
the thermal effect are similar in magnitudes. With further increase of separation an implicit 
thermal effect in the Casimir free energy and pressure becomes larger than an explicit one. 
We have also confirmed the previously known result 18|, [22] that the thermal effect in 
the Casimir interaction between two graphene sheets contributes more than 80% of the 
Casimir free energy and pressure at relatively short separations exceeding 100 and 150 nm, 
respectively, and that the major contribution to both the free energy and pressure is given 
by the electromagnetic held with transverse magnetic polarization. 

In future it would be interesting to extend our asymptotic approach to real frequency 
axis and obtain simple analytic expressions for the polarization tensor, conductivity, density- 
density correlation functions and dielectric permittivity of graphene at nonzero temperature 
for numerous applications other than the Casimir effect. As an example, one could mention 
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calculation of the reflectivity properties of graphene over a wide range of real frequencies. 
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TABLE I: Different constituent parts of the relative thermal correction to the Casimir free energy 
(column 4) arising (column 2) from an explicit dependence of the polarization tensor on temperature 
and (column 3) from an implicit temperature dependence through a summation over the Matsubara 
frequencies are shown at T = 300 K at different separations between two graphene sheets (column 


1). See text for further discussion. 


a (nm) 

l- r M-i r<0) 1 

m 

l^ (0) 1-1-75)1 

m 

175-1-74)1 

\p\ 

5 

0.0312 

0.0082 

0.0394 

10 

0.0922 

0.0347 

0.1269 

20 

0.2013 

0.1068 

0.3081 

30 

0.2730 

0.1734 

0.4464 

50 

0.3492 

0.2680 

0.6172 

100 

0.4115 

0.3807 

0.7922 

150 

0.4303 

0.4289 

0.8592 

200 

0.4387 

0.4552 

0.8939 

250 

0.4432 

0.4717 

0.9149 

500 

0.4511 

0.5064 

0.9575 

1000 

0.4543 

0.5245 

0.9788 


TABLE II: Different constituent parts of the relative thermal correction to the Casimir pressure 
(column 4) arising (column 2) from an explicit dependence of the polarization tensor on temperature 
and (column 3) from an implicit temperature dependence through a summation over the Matsubara 
frequencies are shown at T = 300 K at different separations between two graphene sheets (column 


1). See text for further discussion. 


a (nm) 

|P|-| P (°)| 

\P\ 

/ ,(0> -i/'ol 
l-P| 

\P\-\Po\ 

\P\ 

5 

0.0144 

0.0019 

0.0163 

10 

0.0506 

0.0123 

0.0619 

20 

0.1365 

0.0493 

0.1858 

30 

0.2086 

0.0967 

0.3053 

50 

0.3014 

0.1822 

0.4836 

100 

0.3928 

0.3096 

0.7024 

150 

0.4234 

0.3716 

0.7950 

200 

0.4374 

0.4072 

0.8446 

250 

0.4452 

0.4300 

0.8752 

500 

0.4584 

0.4790 

0.9374 

1000 

0.4637 

0.5051 

0.9688 
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FIG. 1: (Color online) The relative errors of the approximate methods for calculation of the Casimir 
(a) free energy and (b) pressure between two graphene sheets using (line 1) our asymptotic approach 
taking into account an explicit temperature dependence of the polarization tensor in all Matsubara 
terms and (line 2) discarding this dependence in terms with l > 1 are shown as functions of 
separation. 
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FIG. 2: (Color online) The normalized Casimir free energy per unit area of two graphene sheets 
at T = 300 K calculated (the upper line) exactly and (the intermediate line) taking into account 
only an implicit temperature dependence through a summation over the Matsubara frequencies are 
shown as functions of separation. The normalized Casimir energy at zero temperature is presented 
as a function of separation by the lowest line. 
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FIG. 3: (Color online) The normalized Casimir pressure between two graphene sheets at T = 
300 K calculated (the upper line) exactly and (the intermediate line) taking into account only an 
implicit temperature dependence through a summation over the Matsubara frequencies are shown 
as functions of separation. The normalized Casimir pressure at zero temperature is presented as a 
function of separation by the lowest line. 
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